High— order correlations of rich galaxy clusters 



Alberto Cappi 1 and Sophie Maurogordato^ 



1 Osservatorio Astronomico di Bologna, via Zamboni 33, 1-40126 Bologna, Italy 

2 CNRS; Observatoire de Paris-Meudon, LAM, 5 Place J.Janssen, F-92195 Meudon, 
<^ France 



00 



E-MAIL: 
A. CAPPI: cappi@astbo3.bo.astro.it 



m 

00 

o 
o 

O S. MAUROGORDATO: maurogordato@melamc.obspm.circe.fr 

I 

o 



Accepted for publication in the Astrophysical Journal 
July 15, 1994 



Abstract 

We analyse the two-dimensional all-sky distribution of rich Abell and ACO galaxy 
clusters by using counts in cells and measuring the high-order area-averaged angular 
correlation functions. Confirming previous results, we find a well defined hierarchical 
relation between the two and three-point correlation functions, remarkably constant 
with scale. In the angular range 2° < 6 < 4°, the southern sample, limited at bu < 
—40° and including both Abell and ACO clusters, shows a remarkable hierarchical 
behavior up to the sixth order, while northern Abell clusters give positive correlations 
in the same range only up to the fourth order. 

The inferred deprojected values of the 3-D coefficients Sj, where Sj = ^j/^2~ > are 
similar to those measured for the galaxy distribution, and consistent with theoretical 
predictions. These results are confirmed to the 4th order by our analysis of a 3-D sam- 
ple of Abell and ACO clusters. Assuming that selection effects and / or the absence 
of a cluster fair sample are the reason of the difference between the two galactic hemi- 
spheres, and between Abell and ACO clusters, our results indicate that the statistical 
properties of the cluster distribution originate from the underlying galaxy distribution 
and show that the biasing between clusters and galaxies is non-linear. 

Subject headings: Galaxies: clustering — large-scale structure of the universe 



1 Introduction 

The study of the large-scale structure in the Universe is a central issue in cosmology. 
A major challenge for theories of galaxy formation is to explain the observed properties 
of the galaxy distribution. In this field, a statistical approach is natural and fruitful 
(see Peebles 1980); this means, of course, that large, well defined samples are needed. 

Rich clusters represent the highest peaks of the galaxy distribution, and they are 
potentially a powerful tool for the analysis of the large-scale structure. Their use as 
tracers of the luminous matter distribution has been possible after the publication of 
the Abell catalog of clusters (Abell 1958), whose extension by Abell, Corwin & Olowin 
(1989) gives a whole-sky coverage except at low galactic latitude. 

The spatial two-point correlation function £ C c(?") of Abell clusters was measured 
by Bahcall & Soneira (1983) and Klypin & Kopylov (1983), after a first estimate of 
the angular two-point correlation function o> cc (#) by Hauser & Peebles (1973): they 
found that it can be fitted with a power-law comparable to that of galaxies with slope 
7 ~ —1.8, but with a much larger correlation length. In order to explain the different 
amplitudes of the galaxy and cluster autocorrelations, Kaiser (1984) and Bardeen et 
al. (1986) introduced the concept of bias: the selection of rare peaks above some fixed 
threshold in a smoothed Gaussian density field amplifies the correlation function of 
such peaks relatively to the underlying distribution. Nevertheless, this mechanism is 
not sufficient to reconcile the standard CDM scenario with the statistical properties of 
clusters (Bahcall & Cen 1992; Olivier et al. 1993). A better description of a distribu- 
tion requires the knowledge of the higher-order correlation functions. The three-point 
correlation function £ 3 is particularly interesting: in the biased galaxy formation sce- 
nario, for such dense systems as clusters, its decomposition should give a term of order 
3 (Politzer & Wise 1984) which is not observed. The 2-D analysis of Jing & Zhang 
(1989) and Toth, Hollosi & Szalay (1989), and the 3-D analysis of Jing & Valdarnini 
(1991), have shown that the following hierarchical relation holds for Abell and ACO 



(Abell et al. 1989) clusters: 

6 = Q[6(r„)6(r 6 ) + 6(r 6 )6(r c ) + 6(r„)6(r c )] (1) 

with 0.6 < Q < 1.1. This is the same relation previously found for galaxies (Groth 
& Peebles 1977; Davis & Peebles 1977); its generalization at higher orders gives the 
so-called factorized hierarchical models (see Fry, 1984a & b; Schaeffer 1984), where the 
correlation functions at all orders can be written in the following way: 

Um-,rj) = Y,Q { j a) Y,n N -%(r> ] ) (2) 

a ij 

Each term of equation (2) is associated to a graph, points r j correspond to nodes 
and factors ^(fij) correspond to edges between nodes i and j. Such a hierarchy has a 
tree structure, and each distinct tree shape (or topologically distinct graph) is denoted 
by cc, while the sum over ij denotes a sum over relabellings of a given tree a. This 
relation implies an equivalent one in 2-D, where the coefficients Q" are changed into 
the projected qj: 

^....^^g^n^fy (3) 

a ij 

In 2-D, the relation was shown to hold on galaxy catalogs up to J = 4 (Fry & 
Peebles 1978; Sharp et al. 1984). Szapudi, Szalay & Boschan (1992) found that galaxies 
in the Lick catalog continue to follow the hierarchy up to J = 8, and confirmed the 
validity of the relation (1) for clusters selected in the same catalog (see also Borgani, 
Jing & Plionis 1992). 

Hierarchical models represent a subset of the so-called scale-invariant (SI) models 
(Balian & Schaeffer 1989), defined by the scaling property: 

6(Ar 1 ,..,Ar J ) = A-( J - 1 )^ J (r-i,..,r J ) (4) 

This property is verified by a solution of the BBGKY equations for Q < 1. 



The void probability function (VPF), which depends on all the J-point correlation 
functions (White 1979), allowed to confirm the scale-invariance for galaxies (Bouchet 
& Lachieze-Rey 1986; Maurogordato & Lachieze-Rey 1987) and clusters (Jing 1990; 
Cappi, Maurogordato, & Lachieze-Rey 1991). However, the VPF gives global in- 
formation. In a comparison with theoretical predictions, derived from perturbative 
calculations or N-body simulations, we are interested also in a direct knowledge of 
each high-order moment. In SI models, the volume-averaged correlation functions, 
or connected moments, £j = V~ J J v £jd 3 ri...d 3 rj, (we will assume that the average is 
made over a sphere of radius R) can be written as: 

6 = Sjtf- 1 (5) 

where the Sj are independent of R, and can be expressed as an integral over the unit 
volume: Sj = J v d 3 pi, ...,d 3 pj^j(r pi, ...r pj). In the special case of the hierarchical 
models, and assuming that the Qj jCt are independent of cc, then: 

6(ri,...,P J ) = Q J EE nJV_1 6(ra 6 ) (6) 

a ab 

This means that the Sj and the Qj coefficients can be linked by the following rela- 
tion: Sj = FjJ j ~ 2 Qj, where the coefficients Fj depend on the shape of the integration 
volume and on the slope of the two-point correlation function. At order 3, 

s 3 = 3g 3 ^| (7) 

where J 2 = 72/(3 - 7 )(4 - 7 )(6 - 7 )2 7 and J 3 = (3/4tt) 3 / x d 3 x 1 d 3 x 2 d 3 x 3 (\\x 1 - 
^21111^2 — ^H) 7 (Peebles 1980). Recently Gaztahaga & Yokohama (1993) have pub- 
lished the analytical solution of the J 3 integral assuming a power-law two-point cor- 
relation function, and at higher orders numerical estimates are available (Gaztahaga 
1994). Bouchet et al. (1993) have shown that IRAS galaxies follow the hierarchical 
relation (5) up to order 5, and similarly Gaztahaga (1992) has found that also CfA and 
SSRS galaxies follow the above relation up to J = 4. 



In this paper we present estimates of the area-averaged correlation functions for 
the 2-D cluster distribution and we check them using also a 3-D sample, with the aim 
to test directly the validity of the scale-invariance for galaxy clusters, and to compare 
the high-order moments with those measured for the galaxy distribution and expected 
from the theory. The method is presented in section 2; section 3 is devoted to the 
definition of the samples and the presentation of the main results; in section 4 the 
implications of these results are discussed. Our main conclusions are summarized in 
section 5. 

2 Counts in cells and moments 

We calculate the moments of the cluster distribution from counts in cells (see Bouchet, 
Davis & Strauss 1992, Gaztahaga 1992, and Bouchet et al. 1993). Here we describe 
the procedure for the analysis of a 2-D catalog. Given a sample of objects with known 
angular coordinates, we place randomly a fixed number of circles with radius 8 and 
area 2ir[l — cos(#)]. A number of circles N ran = 10 6 was found sufficient to avoid 
undersampling at small scales and to have negligible fluctuations. Circles which cross 
the sample borders are not counted. We measure the count probabilities P(N) of 
finding N objects within a radius 8, by counting the number of circles containing N 
objects and dividing it by the total number of circles. This procedure is repeated at 
different scales 8. 

The same method can be applied to the analysis of 3-D volume-limited samples; 
instead of circular areas, spherical volumes of radius R are used. From the P(N) we 
can estimate the centered moments of order J (see Peebles 1980): 

where N = J2 NP(N), i.e. at each scale we compute directly the effective N from 
counts. This means that the effective N may slightly differ from the sample N = nV, 



where n is the mean density of the sample. We choosed to use the effective N for 

self-consistency; anyway we verified that the two methods of computing N give similar 

results. 

The area (or volume) -averaged correlation functions (the reduced moments) are 

finally given by the following relations, where the Poisson shot noise is subtracted from 

each moment: 

1 

^2 = 1*2 ~ ft (9) 

-3^3-3^ + J, (10) 

C^ 3 I 11 ^ 2 9 2 6 

"4 = iia ~ 6— + 11^ - 3/x 2 - -^ (11) 

tf 6 = /* 6 - 10^ + 35^ - 50^ + 30^ - 1(W 3 + ^ (12) 

^e = ^,-15^+85^-225^+274^-255^+30^-10^ + 150^-15^-^ 

(13) 

3 Analysis of the samples 
3.1 Definition of the 2— D samples 

We analysed both 2-D and 3-D cluster samples. The advantage of using 2-D samples 
is the larger number of objects, which allows a better estimate of the higher moments; 
but in this way one can only indirectly obtain the corresponding 3-D values, and for 
this reason the analysis of smaller 3-D samples is a complementary and useful check. 
The 2-D sample includes 1674 clusters (NS45D61), 774 Abell clusters in the north 
(N45D61) and 900 Abell + ACO clusters in the south (S45D61) -hereafter for Abell 
clusters we will mean the clusters at 8 > —27° selected by Abell (1958), and for 
ACO clusters the southern clusters selected by Abell, Corwin, Olowin (1989)- with 
I bu |> 45°, richness class R > 1 and distance class D < 6, covering an area of 3.68 
steradians (see table 1). The high cut in galactic latitude allows to minimize the effect 



of galactic extinction (we analysed also a sample limited at | bu |> 40°, without finding 
significantly different results). 

For what concerns projection effects, Olivier et al. (1990) found that projection 
contamination affects significantly o> cc (#) only at scales below 1°. Peculiar motions 
of clusters are small relatively to their Hubble velocities (Huchra et al., 1990), and 
distortions in redshift space should be negligible, as well as their effect on the esti- 
mates of Sj. Finally, the southern samples are a mixture of Abell and ACO clusters; 
while having been selected according to the same criteria, they show some differences 
(see Scaramella et al. 1991) which nonetheless might be partly real (e.g. Cappi & 
Maurogordato 1992). 

Errors were computed by the bootstrap technique (see Ling, Frenk & Barrow 1986). 
We generated 50 catalogs from each sample, then we estimated errors in two ways. 
For each catalog, we found the values of sj and -jj from a linear fit of the relation 
log a) j = y4 + 7jloga>2 in the range 2° < 8 < 4°; then we estimated the respective rms. 
We also fixed -jj = J — 2, and we estimated the rms of sj at different angular scales. 

3.2 Analysis of the 2— D samples 

Figures la-d show u)j, J = 3, ..,6, as a function of u) 2 - Looking first at the total sample 
(black hexagons), a very good hierarchical behavior appears, especially in the range 
2° < 6 < 4°. It is noteworthy to stress that the straight lines in the figures are not fits 
of the points: they simply represent the relation ujj = J j ~ 2 uj 2 ~ ? which differs from 
the hierarchical one by a factor of the order of unity (e.g. Gaztahaga 1994). Table 2 
lists the average values of sj calculated assuming the hierarchical relation (5) between 
2° and 4°, and also the values of -jj and sj computed from a least-square fit to the 
relation log(u)j) = log(sj) + / yj\og(u>2)i leaving the slope -jj as a free parameter, with 
the corresponding bootstrap errors. Notice that for the fit we kept the range fixed 
(2° < 6 < 4°) at all orders. 



Figures 2 and 3 show respectively uij and s j as a function of scale. The correlations 
are consistent with power-laws systematically steepening at higher orders. s 3 , s 4 , 
5 5 and s 6 are nearly constant between 2° and 4°, corresponding roughly to a spatial 
range of 20-40 h~ x Mpc. The analysis is limited at small scales by shot noise, and 
at large scales by the small number of independent areas, with an increasing error in 
the estimate of P(N) at large N, and by the loss of the P(N) exponential tail, with a 
systematic underestimate of the higher moments. Moreover, at large scales systematic 
discrepancies appear between Abell and ACO clusters (see Olivier et al. 1990). 

Up to now, we have discussed the total sample, stressing its hierarchical behavior. 
But the situation appears more complex when looking separately at the northern and 
southern samples, because a systematic discrepancy is apparent. In figures la-d, it is 
possible to see that most of the hierarchical behavior is due to the southern sample. 
Northern Abell clusters follow the hierarchical behavior at order 3 (with a lower 53 
than southern clusters), but they deviate from the expected behavior at the 4th order, 
and their u> 5 becomes negative very rapidly, while the S45D61 sample shows sj nearly 
constant in the range 2° < 8 < 4° up to the 6th order. These effects are apparent 
in figure 4, which shows sj as a function of 8 for the northern and southern samples 
separately. 

That a difference could be detected is not surprising. For example, 1^2(8) of southern 
clusters [bu < —45°) is lower than the northern one at large angles (as found for the 
standard 2-point correlations by Olivier et al. 1990), while in the range 2°-4° the 
values of s j for southern clusters are systematically higher than the northern ones. 
These effects become more important not only at larger 8, but also at higher orders. 

Such systematic differences constitute the main limit of any inference from the Abell 
and ACO cluster distribution. 

In principle, according to bootstrap errors, the difference between the northern 
and southern sj at high orders is significant, but errors derived with the bootstrap 
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technique, while more realistic than Poissonian ones, might be an underestimate. 

Moreover, it is probable that we do not dispose yet of fair samples. The presence 
of some large structures might change drastically the estimated count probabilities, 
especially at the tail of the distribution. This might explain the relative insensitivity 
of the two-point correlation function or the VPF to these differences, which become 
increasingly large at higher orders. 

Also selection effects in the Abell and ACO cluster catalog might cause some dif- 
ferences, and the ACO catalog might be deeper than the Abell one (see e.g. Plionis et 
al. 1992). We checked this possibility analysing two subsamples, the first made only of 
Abell clusters at southern galactic latitude [bu < —40°, 8 > —17°), the second made 
only of ACO clusters {bu < —40°, 8 < —27°). The results are quite noisy, especially 
at the higher orders, nevertheless they show a systematic difference between ACO and 
Abell clusters, and only ACO clusters show clear signs of a hierarchical behavior at 
orders > 4. 

We suspect that a combination of the above discussed effects is at the origin of this 
discrepancy. Nevertheless, it is interesting to note that the combination of these sub- 
samples gives a stronger hierarchical behavior; this fact gives us some confidence that 
cluster correlations are indeed hierarchical, but that in present catalogs this hierarchy 
appears clearly only at the lowest orders, where the signal-to-noise is sufficiently high. 

3.3 Deprojection of sj 

We have obtained the projected sj, and the next step is the estimate of the 3-D Sj. 

The values found for the total sample NS45D61 -assuming that they represent our 
best estimates of the correlation hierarchy (5) at scales between 2° and 4° degrees, cor- 
responding approximately to the range 20 - 40 h~ x Mpc- are the following: s 3 ~ 3.13, 
5 4 ~ 16, s 5 ~ 110, s 6 ~ 745. The derivation of Qj from sj requires two steps. Assum- 
ing that the the spatial correlation functions £j follow the hierarchical relation (5), a 



2-D projection gives the corresponding relation for the angular correlation functions 
u)j, where the 2-D qj differ from the 3-D Qj by a multiplicative factor rj, qj = rjQj: 



a 3 



gf^n^^) = P 3 Q ( f ] rrldr^in) f°° r 2 2 dr 2 $(r 2 ) 
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In equation (14) $(r) is the selection function of the sample, a and p are respec- 
tively the 2-D and 3-D densities. This equation can be solved analytically for small 
separations, assuming a power-law for the two-point correlation function. 

Toth et al. (1989), assuming a step-like form for the selection function and a power- 
law form for the two-point correlation function, have calculated analytically the ratio 
between Q 3 and q 3 , finding Q 3 = 5 3 /1.040. Bernardeau (1994, private communication) 
has extended the calculation up to the 6th order. At this point, ujj can be computed 
from the area-averaged uij. At each order, this implies the computation of an fj 
integral and a sum over all possible vertices, leading to sj = fjJ J ~ 2 qj. Numerical 
estimates of the fj integrals (Bernardeau 1994, private communication) show that 
the corrective factors are around unity, and not very sensitive to the geometry or to 
variations of the £ 2 slope. In the same way, the 3-D Sj and Qj are related through 
Sj = FjJ j ~ 2 Qj, where the Fj are integrations over the 3-D box. Fj have been 
calculated analytically at order 3 (Gatzanaga and Yokohama 1993) and numerically 
at higher orders (Gatzanaga 1994) in the case of a power-law slope with 7 = —1.8 for 
the correlation function. In this way we can obtain a first rough estimate of the values 
of Sj deduced from the 2-D analysis. Table 3 gives the main quantities involved in 
this computation, as well as the values of Qj and Sj. The estimated values of Qj 
for clusters are Q 3 = 1.0, Q4 = 0.9, Q$ = 0.7, Qe = 0.4. Detailed calculations of 
these factors, as well as a more complete analysis using a selection function computed 
from the data will be presented in another paper (Bernardeau, Cappi, & Maurogordato 
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1994). 

3.4 Analysis of the 3-D sample 

We checked our results analysing a nearly complete spatial sample limited at | bu |> 40° 
and z < 0.08, with 281 clusters (120 Abell clusters in the north and 161 Abell + ACO 
clusters in the south, including richness class 0). It is an extension of the sample defined 
by Cappi & Maurogordato (1992), without the limit on the apparent magnitude of the 
tenth cluster member (m < 16.4) derived from the Postman et al. sample (1992). 
Comoving distances in Mpc are calculated from redshifts using the Mattig formula: 



c 1 

H~o7o' 



D = —-[1 - q + q z + (q - l)j2q z + 1]/(1 + z) (15) 



where we fixed H = 100 and q = 0.5 (the only choice which allows using euclidean 
geometry for separations and volumes). 

The 3-D sample is too small to give significant results for the moments higher than 
4, and we cannot test if the discrepancy found in the analysis of the 2-D samples is 
reproduced here. Considering the total sample (both hemispheres), we obtain ^3 ~ 2.5 
between 22 and 38 h _1 Mpc, and S4 ~ 12 between 30 and 35 h _1 Mpc (fig. 5). 

Here again, we find a systematic difference between the northern and southern 
samples; in particular, we cannot measure a positive £4 for the northern sample, except 
for a small range around 8 ~ 2 degrees. It may be interesting to note that the northern 
and southern samples are different also in their percolation properties, and that this can 
be explained by the presence of the large Horologium-Reticulum supercluster (Cappi 
& Maurogordato 1992). Therefore it is possible that the systematic difference, which 
becomes significant at high orders (see also the multifractal analysis of Borgani et al., 
1994), is partly due to a real difference in the cluster distribution. The 2-D and the 
3-D analysis give consistent values of Sj and comparable scales where the hierarchical 
behavior is observed. This is noteworthy, because the 3-D sample is limited at z < 0.08 
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and includes richness clusters, while the 2-D samples are mainly composed by clusters 
of distance class 5 and 6 and richness class R > 1. We can conclude that clusters follow 
a hierarchical behavior up to the 4th order, and -from the analysis of the 2-D samples- 
that southern clusters follow a hierarchical behaviour up to the 6th order. 

4 Discussion 

Recovering the 3-D coefficients Qj and Sj allows a comparison between our results 
and those obtained for different galaxy samples, and also with theoretical predictions. 
From our 2-D analysis we find Q 3 ~ 0.98, comparable to the values found with other 
methods. Gaztahaga (1994) estimated Sj -up to the 9th order- from the galaxy 
angular distribution in the APM Survey. At scales between 7 and 30 h _1 Mpc, he finds 

5 3 = 3.1 ± 0.14, S 4 = 22.0 ± 2.5, S 5 = 146 ± 47, and S 6 = 1514 ± 337. A comparison of 
these values with the corresponding cluster values reported in table 3 shows that the 
galaxy and cluster ^3 values are remarkably similar, while their 5s are consistent at 
the 1-c level, their S4 values differ at the 2-cr level, and only Se values are different 
at the 3-cr level. While uncertainties (depending also on the deprojection factors) do 
not allow a detailed comparison, these results suggest that the APM Galaxy Survey 
and rich clusters are sampling the same underlying large-scale distribution. In the 
hierarchical tree models studied by Bernardeau and Schaeffer (1992), for high density 
objects such as galaxies or clusters, the values of all Qj(x) converge to 1 for large x 
-and even as soon as x > 1- (where x = N/N c is the scale parameter). The values 
of Q j for clusters are in good agreement with these models, except for Q 6 , which 
however has the highest error, as clearly indicated by the large bootstrap errors on 
s 6 and 7 6 . The values of Sj are also similar to the predictions derived in the frame 
of perturbation theory; assuming a self similar power spectrum and a spectral index 
n = —1, one gets ^3 = 2.86 (Juskiewicz, Bouchet, Colombi 1993; Bernardeau 1994) and 

54 = 23 (Bernardeau 1994; Catelan & Moscardini 1994); from numerical simulations 
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with scale-free Gaussian initial conditions, Lucchin et al. (1994) find (again for n = — 1) 
^3 ~ 3.1 — 3.3, 54 ~ 16 — 18 and S 5 ~ 135 — 147, while in their N-body simulations of 
a Cold+Hot Dark Matter universe Bonometto et al. (1994) find S 3 ~ 2.5 and S^ ~ 7.5 
for galaxies. 

Following Fry and Gazahaga (1993), the smoothed biased (galaxy) density field can 
be expressed as a Taylor series of the matter density field: 



oo 



Sb = £ t/m (16) 

fc=o K - 

In the linear approximation, we have Sb ~ Wm and ^ 2 fc ~ b 2 £ 2 M- It follows that the 
"normalized" skewness S3 should be inversely proportional to the bias: 



S 3b =< (bS M f > I < (bS M ) 2 > 2 = S 3M /b (17) 

If we take into account higher orders in 8m in the development of Sb, the relation 
between £j and £ 2 is affected. For instance, taking only the leading term in ( 2 M, we 
have for S3: 

S 3b = b- 1 (S 3M + Sb 2 /b) (18) 

while new terms appear at higher J (Fry & Gaztahaga 1993; Bernardeau 1994). 

For J = 2, the linear bias expression for £ 2 coincides with the leading term in the 
non-linear bias expression, and we have £ 2 fc ~ b 2 £ 2 M in both cases. From this relation, 
it is possible to measure for example the relative bias between clusters and galaxies 
(see Peacock k Dodds 1994). 

For J = 3, the linear biasing -equation (17)- and the non-linear biasing -equation 
(18)- give two different relations, and this allows to discriminate between the two cases. 

We have shown that galaxies and clusters have nearly the same S3: as a consequence, 
a linear biasing is ruled out. 
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Moreover, the fact that the values of Q j are around unity for both galaxies and 
clusters is consistent with models where a natural bias comes out from gravitational 
interaction of the matter field. 

5 Conclusions 

We have calculated the high-order connected moments (up to order 6) of the rich 
cluster distribution using counts in cells. Our main results are based on the analysis 
of the 2-D distribution, and on a smaller 3-D sample. 

We have shown that the all-sky 2-D and 3-D cluster distribution is consistent with 
the correlation hierarchy (5) up to J = 4, and the 2-D southern cluster distribution 
up to J = 6, at scales between 2° and 4° degrees, corresponding approximately to the 
range 20 - 40 h' 1 Mpc. 

The inferred values for Sj are comparable to those found for galaxy samples; in 
particular, the similar values of ^3 for galaxies and clusters imply a non-linear bias. 

Anyway, subsamples show systematic differences. Northern Abell clusters show a 
clear hierarchical behavior only at order 3. These differences can be due to a com- 
bination of selection effects, smallness of samples and differences in the real cluster 
distribution. The fact that the combined, largest sample shows a well defined hierar- 
chical behavior even at the 6th order, suggests that cluster correlations are hierarchical, 
but only larger samples will allow to confirm this result. 

These statistical properties may give stronger constraints to the various models of 
galaxy formation even if the accurate definition of clusters in N-body simulations on 
one hand, and the construction of unbiased, fair cluster samples on the other hand, 
represent problems yet to be solved. 
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Figure captions 

Fig. 1 a— d Area-averaged angular correlations ujj (J = 3,4,5,6) as a function of the two- 
point correlation u) 2 - for the total sample NS45D61 (black hexagons), the north- 
ern sample N45D61 (open squares) and the southern sample S45D61 (open hexagons). 
Solid lines are power-laws with the following parameters: a) s 3 = 3, 73 = 2; b) 
54 = 16, 74 = 3; c) 5 5 = 125, 75 = 4; d) s 6 = 1296, 7 6 = 5. 

Fig. 2 Area-averaged angular correlation functions u)j, J = 2,3,4,5,6, as a function 
of scale 8, for the total sample NS45D61. Lines represent least-squares fit of 
the points, log(a>j) = Aj + Bjlog(8), within the range 2° < 8 < 4°. Best-fit 
parameters are the following: A 2 = 0.033, B 2 = -0.95; A 3 = 0.538, B 3 = -1.84; 
A 4 = 1.275, £ 4 = -2.79; A 5 = 2.183, B 5 = -3.81; A 6 = 3.093, B 6 = -4.86. 

Fig. 3 Values of sj = ujjIuj{~ x (J = 3,4,5,6) for the total sample NS45D61. Bootstrap 
errors are shown in the range 1° < 8 < 5°, with step equal to 1°. 

Fig. 4 Values of sj = u>j/u> 2 ~ 1 for the S45D61 (open points) and N45D61 (black points) 
samples. J = 3: triangles, J = 4: squares, J = 5: pentagons, J = 6: hexagons. 
Bootstrap errors are shown in the range 1° < 8 < 5°, with step equal to 1°. 

Fig. 5 Values of Sj = ^j/^ 1 {J = 3,4), resulting from the analysis of the 3-D cluster 
sample. Bootstrap errors are shown in the range 20-40 h~ x Mpc, with step equal 
to 5 h' 1 Mpc. 
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Table 1: Definition of the 2-D samples. 





N45D61 


S45D61 


bn 


±45° 


-45° 


N c 


774 


900 


-^Wnm 


1 


1 


Umax 


6 


6 


Z"max 


— 


— 



Table 2: Average values and rms of sj between 2° and 4°, and values of sj and -jj 
estimated from a linear fit of the relation log u)j-log u) 2 in the same range, with the 
corresponding bootstrap errors. 



J 3 


4 


5 


6 


sj (mean and rms) 3.13 ± 0.06 
sj (fit) 2.981°;^ 
7 j 1.95 ±0.09 


16 ±1 

i5.2±f : | 

2.96 ±0.34 


110 ±12 
4.04 ±1.37 


745 ± 123 

855t^ 6 
5.16 ±3.84 



Table 3: Factors used for deprojecting the 2-D estimates of sj. 



J 


3 


4 


5 


6 


sj 


3.13 


16 


110 


745 


fj 


1.02 


1.03 


1.08 


1.15 


Fj 


1.03 


1.10 


1.15 


1.25 


rj 


1.04 


1.11 


1.22 


1.36 


Qj 


0.98 


0.87 


0.67 


0.37 


Sj 


3.04 


15.4 


96 


595 
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Figure 1a 
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Figure 1c 
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Figure 1d 
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